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Abstract 


Time evolution of plane, cosmic-ray modified shocks has been simulated numerically 
for the case with parallel magnetic fields. Computations were done in a ’’three-fluid 
dynamical model incorporating cosmic-ray and Alfven wave energy transport equations. 
Nonlinear feedback from the cosmic-rays and Alfven waves is included in the equation ot 
motion for the underlying plasma, as is the finite propagation speed and energy dissipation 
of the Alfven waves. Exploratory results confirm earlier, steady state analyses that found 
these Alfven transport effects to be potentially important when the upstream Alfven speed 
and gas sound speeds are comparable. As noted earlier Alfven transport effects tend to 
reduce the transfer of energy through a shock from gas to energetic particles. These studies 
show as well that the time scale for modification of the shock is altered in nonlinear ways. It 
is clear, however, that the consequences of Alfven transport are strongly model dependent 
and that both advection of cosmic-rays by the waves and dissipation of wave energy in 
the plasma will be important to model correctly when quantitative results are needed. 
Comparison is made between simulations based on a constant diffusion coefficient and 
more realistic diffusion models allowing the diffusion coefficient to vary in response to 
changes in Alfven wave intensity. No really substantive differences were found between 

them. 
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1. Introduction 


Diffusive transport of cosmic-rays in shocks is mediated through resonant scattering 
by Alfven waves. The dominant interactions probably involve Alfven waves stimulated to 
propagate in the direction of cosmic-rays streaming away from the shock itself. Although 
the existence of such waves is clearly implicit in all discussions of diffusive shock acceler- 
ation of cosmic-rays, most treatments of this problem have neglected the motion of the 
Alfven waves through the background plasma and the energy exchange between the Alfven 
waves and the plasma; what we may term ’ Alfven transport effects. Several authors have 
pointed out, however, that Alfven transport effects can place important controls on particle 
acceleration and shock dynamics when the Alfven speed is significant compared to the up- 
stream gas sound speed (e.g., Holman, Ionson k Scott 1979; Achterberg 1982, Volk, Drury 
and McKenzie 1984). To date only steady shocks have been discussed in this context. As 
a first step to examine these issues in unsteady cosmic-ray shocks I have extended the 
time dependent two-fluid numerical methods for diffusive cosmic-ray transport employed 
in Jones k Kang (1990) to include Alfven transport in parallel MHD, sonic-mode shocks. 
More nearly complete Alfven transport formalisms based on the diffusion- ad vection equa- 
tion for the cosmic-rays and general MHD are under development. In the present case the 
governing equations for the gas are: 


— = — pV ■ it. (1-1) 

dt 1 

( IlL = -Iv(P v + Pc + Pw), (1-2) 

dt p 

and 1 £ 

%■ = --V ■ {(P, + Pc + Pw)Z} + ~(Pc + Pw)$ • S + (1-3) 

dt p P P 

where u is the velocity of the gas. while P ,j , P c — (lc — 1 )E C an d P-w — E w /2 are the 
pressures of the gas, the cosmic-rays and the waves, respectively, d/dt is the Lagrangian 
time derivative including the gas motion, e is the sum of gas thermal and kinetic energy 
per unit mass. The term L represents nonadiabatic heating from dissipation of the Alfven 
waves. Equations [1.1 - 1.3] differ from those in the usual two-fluid model for cosmic-raj 
transport only through the addition of the wave pressure, P w . 


The cosmic-ray energy density is followed through an energy conservation equation 
derived from the diffusion-advection equation given by Skilling (1975). It corrects the 
standard two-fluid expression (e.g., Drury and Volk 1981) for the advected cosmic-ray 
energv flux and the work done by the cosmic-rays on the gas for the fact that the cosmic- 
rays are advected with the waves rather than the gas (McKenzie and Volk 1982, Achterberg 
1982): 


dE c 
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where U) is an energy weighted, mean diffusion coefficient and u w is the center-of- 
momentum velocity for the resonant Alfven waves. Equation [1.4] is valid when this 
Alfven center-of-momentum frame is independent of wavenumber, so that its motion is 
”fluid-like”. If the wave propagation is unidirectional, then u w ~ v A in the direction o 
wave propagation. In a more general case, u w would lie between ±v A . For the latter case 
one should properly include an energy gain term in equation [1.4] (and an analogous loss 
term in the energy equation for the Alfven waves) to account for second order Fermi ac- 
celeration; i.e, momentum diffusion. However, within the shock structures energy transfer 
due to momentum diffusion will generally be smaller than the spatial diffusion term y a 
factor less than ~ ( v A /u s ) 2 , where u s is the shock speed relative to the unshocked gas. bo 
long as we limit ourselves to superalfvenic shocks, this should not have a major impact on 
energy transfer, although it might be significant in determining the form of the cosmic-ray 
momentum distribution. 

Two-fluid methods have the limitations that they cannot follow the evolution of the 
cosmic-ray particle distribution self-consistently and that, consequently, one must a prion 
specify closure parameters such as the cosmic-ray adiabatic index, 7c, and the mean 
diffusion coefficient, (k), although they need not be constants. On the other hand, within 
the bounds implied by those caveats two-fluid models accurately describe the dynamics 
of the shock evolution and provide an efficient means to study the approach to dynamical 
equilibrium in such shocks. Our previous calculations comparing two-fluid and diffusion- 
advection calculations have shown two-fluid models to be an effective means to explore 
the general dynamical evolution of cosmic-ray mediated shocks (Kang k and Jones 1991, 
Kang, Jones & Ryu 1992, Jones Kang 1992). 

To complete the formalism an energy equation for the Alfven waves is needed. In the 
WKB approximation and under constraints similar to those applied to equation [1.4] the 
wave energy density obeys 

^ = — i7 u . ■ VPr - V • (u w E w ) - *E U V • u - L. (1.5) 


A full discussion of equations [1.1 - 1.5] with associated references is given in Jones 
(1993) [hereafter J93]. That paper also explores a range of plausible properties for Alfven 
transport as applied to the evolution of plane shocks with magnetic fields parallel to the 
shock normal. The key Alfven properties are; a) the model for the effective wave transport 
velocity, i7 w , as a function of space and time, b) the model for wave energy dissipation 
I c) the ratio of the upstream value of the Alfven speed, v Ao , to the upstream gas sound 
speed. c so (v Ao = c S 0 /y/P), and d) the ratio of the upstream wave energy density to the 

magnetic energy density, a = E W o! E Bo- 

Very briefly, the results discussed in J93 confirm the findings of earlier, steady state 
calculations showing the potential importance of Alfven transport effects (e.g., Volk, Drury 
and McKenzie 1984). Generally, Alfven transport reduces the rate for transfer of energy 
from gas to cosmic-rays and reduces the time-asymptotic post-shock cosmic-ray pressure^ 
For shocks of moderate strength (sonic Mach numbers up to a few tens), for values of 
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3 ~ 1, and for u w ss the time-asymptotic, post-shock cosmic-ray pressure can be 
reduced by more than an order of magnitude by these effects. An example illustrating this 
impact is shown in Figure 1. Details of those models are described in the next section and 
more completely in J93. Thus Alfven transport has the potential to practically remove 
the strong modification of the shock structure by the cosmic-rays. The time necessarv 
for a modified shock to approach dynamical equilibrium is also changed, but in nonlinear 
wavs that sometimes increase and sometimes decrease the equilibrium time compared 
to that found in similar calculations with negligible Alfven transport. It was clear from 
the time dependent simulations that influences both from the advection by waves and 
from wave energy dissipation are important contributoi's to the dynamical evolution of 
cosmic-ray modified shocks. The net influence of Alfven transport depends sensitively 
upon the details of the generation, dissipation and especially advection of Alfven waves in 
the vicinity of shocks. In addition, it was clear that Alfven transport properties on both 
sides of the shock are important whenever one considers unsteady shocks. For example, 
in the event that streaming by cosmic-rays downstream of the shock transition stimulates 
Alfven waves propagating downstream as well, the net effects of Alfven transport can be 
greatly increased relative to models in which the waves always face upstieam. 


2. Evolution of the Mean Diffusion Coefficient 


One effect that was not discussed in detail in J93 was the influence of the spatial 
and temporal evolution of the diffusion coefficient, (k) due to changes in the Alfven wave 
intensity according to equation [1-5]. In theory the diffusion coefficient depends upon the 
Alfven wave energy density, according to 


where Eg = B 2 /{8 tt ), is the total magnetic field energy density. Since the calculations 
described in J93 showed increases in E w by factors ~ 10 across strong shocks, whereas 
Eg remains constant in parallel shocks, equation [2.1] would imply large decreases in 
(k) as material passes through. In the simulations discussed in J93 (k) was taken to 
be a constant in order to isolate other important physical effects. The question would 
be what impact various properties of (k) have on the dynamical evolution of the shocks 
under study. On the one hand, (k) determines the characteristic length scale for the 
shock precursor, xj ~ {k)/u s , and the timescale. t d = x d /u s , for dynamical evolution 
of a shock; i.e., the time necessary to approach an approximate dynamical equilibrium. 
Those scales will be reduced as ( k ) becomes smaller, as will the time for individual cosmic- 
ray particles to be accelerated. On the other hand, simple arguments suggest that the 
value of (k) cannot influence steady state solutions to equations [1.1-1. 5]. Of course, 
on sufficiently long timescales the details of the diffusion coefficient, and especially its 
momentum dependence must influence such issues as the eventual escape of very high 
energy cosmic-ray particles. For that reason and others, as well, steady two-fluid models 
are probably not very realistic as detailed models. Strongly mediated cosmic-ray shocks 
are perhaps never in true equilibrium, anyway, since they seem to be subject to various 
instabilities (e.g., Ryu, Kang k Jones 1993). Thus the specific question more relevant to 
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the present explorations is what influence the properties of (k) in equation [ 2 . 2 ] may have 
on the timescales and on the properties of cosmic-ray mediated shocks as they approach 
dynamical equilibrium. 

To allow simple examination of these questions. I present results from three shock 
simulations (labeled below as models A, B and C ) that differ only in the assumed character 
for the mean diffusion coefficient, (k). The other, common properties are as follows: All 
the models assume a uniform initial (upstream and downstream) cosmic-ray pressure, 
p t) _ p (Q — i 35 x 10 “2, a gas Mach number for the initial shock M = 10 , a ratio of 
upstream sound speed to Alfven speed. 3 = 1. and an upstream ratio of energy density 
in Alfven waves to total magnetic field, a = 0.5. The cosmic-ray adiabatic index is 
assumed to be = 1 . 4 . which is approximately what one would infer for galactic cosmic 
rays if the observed power law momentum distribution extends to nonrelativistic energies. 
That choice is not critical to our results, however. I assumed in each simulation that 
| t 7 u ,| . VA ' with the sign determined by the local cosmic-ray streaming direction; i.e, so 

that u w • VP C < 0. I also assumed a "local equilibrium” model for the wave dissipation 


All the simulations were initiated with a preexisting gas shock discontinuity propagating 
to the right at unit speed in the frame of the grid. That corresponds in the frame of the 
upstream gas to u s = 1.5. Numerical methods were those described in detail in J93. 

For model A the diffusion coefficient was taken to be the constant (k) = 0.35. This 
model is thus identical to the model 3 in J93. except for the value of a. For model B I took 
(„•) = 1 . 4 /p. while for model C (k) = IAE W o/E u . The behavior for (k) in model B follows 
that used by a number of authors before (e.g.. Drury & Falle 1986. Jones and Kang 1990 ). 
It. was originally introduced to avoid the possible development of an acoustic instability in 
strongly cosmic-ray modified shocks, but also provides a crude way to allow for the effects 
indicated by equation [2.1]. Model C is. of course, directly based on equation [2.1]. 

The evolution of all three shocks A - C is illustrated in Figure 2. It is readily apparent 
that the differences among them are. in fact, minor. Because (k) becomes smaller in 
models B and C, the effective length and time scales, xj and are reduced. T us 
the width of the shock precursor is smaller, and the shocks reach dynamical equilibrium 
more quickly. The former effect is best seen in the plot of P c , while the latter is evident 
in all the variables, since it shows equilibrium postshock properties farther behind the 
"current” shock positions. For all times shown, however, all the shocks have actually 
closely approached equilibrium in their immediate postshock properties. 

On the other hand, and perhaps most important to the issues being addressed in the 
study by J 93 , all three models produce very similar time asymptotic values for P g and P c \ 
i.e, the ” efficiency” of particle acceleration and the degree of shock modification by cosmic- 
rays is the same' to within the limits of the numerical scheme (see J93, for a discussion 
of' those limits). Thus, for qualitative studies of this type, involving plane shocks and 
two-fluid models, the detailed behavior of <k) does not seem to be critical. However, for 
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unsteady, nonplanar shocks, such as supernova remnants (e.g., Jones & Kang 1992, Jones 
& Kang 1993) the specific properties of (k) are more significant. 

3. Conclusions 


Alfven transport effects in ’’low 0” plasmas can be important to the transfer of energy 
from the thermal plasma to high energy cosmic-ray particles. Because streaming by the 
cosmic-rays tends to follow the propagation of Alfven waves resonantly interacting with 
them, and because the speed of that propagation in parallel shocks tends to change across 
the shock transition, both the energy gain by a particle and the residence time by particles 
near the shock are modified from calculations which assume that cosmic-rays are advected 
exact lv with the gas. Furthermore, dissipation of wave energy in the gas transfers energy 
from cosmic-rays to the gas and preheats gas entering a shock, tending for both reasons to 
reduce the net efficiency of energy transfer to cosmic-rays in the shocks. In an associated 
paper ( J93) I have discussed in detail these effects in the context of the dynamical evolution 
of plane shock structures, emphasizing the sensitivity to uncertain details about models for 
Alfven transport. Those computations are based on a three-fluid model for the dynamical 
interactions of the gas, the cosmic-rays and the Alfven waves in shocks with magnetic fields 
parallel to the shock normal. That discussion did not consider consequences of using a. 
self-consistent model for the diffusion coefficient which allowed for its dependence upon 
the local intensity of Alfven waves. Rather, in an effort to isolate other, clearly important 
model features for Alfven transport, J93 presented only calculations using constant mean 
diffusion coefficients. In the present paper I have presented a comparison of one such model 
with models that allow for changes in the diffusion coefficient in response to changes in the 
Alfven wave energy density as the waves are advected through a shock. These computations 
are still done within a three-fluid model, but in one case do include nonlinear modifications 
to the waves due to the cosmic-rays and to wave dissipation within the plasma. Although 
there are expected changes in the characteristic diffusion lengths and times associated with 
the shock structure, the more complicated behavior of the diffusion coefficient does not 
appear to alter in any discernible way the time-asymptotic properties of the postshock 
flow. Most important in this regard, it does not seem to alter the expected net efficiency 
of energy transfer to cosmic-rays within a transition. 


Of course, treatment of diffusive shock acceleration with Alfven transport based on 
a more complete formalism using the momentum dependent, diffusion-advection equation 
will be necessarv before one could possibly extend that conclusion to any discussion of the 
properties of the cosmic-ray particle distribution. To that end I have under development an 
effort to simulate the evolution of cosmic-ray mediated shocks using the diffusion-advection 
equation and also including nonlinear magnetohydrodynamical (MHD) gas dynamics (see 
Frank. Jones & Ryu 1993) (this proceedings]. The much greater complexity of those 
treatments, however, will certainly cause us to look back for bench-marks to simpler 
models, such as those discussed in J93 and the present paper. 
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Figure Captions 


Pig, i. — Evolution of two plane shocks differing only in the presence or absence of Alfven 
transport effects. The solid curves represent model A, described in the text and 
also shown in Figure 2. The dotted curves represent the behavior for the same 
model with 3 = 10 4 . so that the Alfven transport, terms in equations [1.1 - 1.4] 
are negligible. Density, gas velocity, gas pressure and cosmic-ray pressure are shown 
at times t = 0.30.60,90. The nominal diffusion times and lengths for these models 
are t d = 0.156 and x d - 0.233. Note, in addition to the large change in P c between 
the two models the substantial change in compression and shock speed. 

pig 2. — Evolution of three shocks using different models for the diffusion coefficient. 
(k). Each simulation assumes ~ c ~ 1-4, 3 = 1, a = 1/2, M = 10, P C o — Pgo ■ The 
solid curve represents the case (k) = 0.35 (model A). The dotted curves correspond 
to (k) = 1.4( E w 0 /E w ) (model C). The dashed curve ( P c only) shows the case ( k ) = 
1.4/p (model B). The asymptotic shock properties are the same to within the limits 
of numerical accuracy at finite time and resolution. Differences in P c within the 
precursors are due to differences m the effective diffusion length, x d . 
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